Mobility of Discrete Solitons in Quadratically Nonlinear Media 
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We study the mobility of solitons in second-harmonic-generating lattices. Contrary to what is 
known for their cubic counterparts, discrete quadratic solitons are mobile not only in the one- 
dimensional (ID) setting, but also in two dimensions (2D). We identify parametric regions where 
an initial kick applied to a soliton leads to three possible outcomes, namely, staying put, persistent 
motion, or destruction. For the 2D lattice, it is found that, for the solitary waves, the direction 
along which they can sustain the largest kick and can attain the largest speed is the diagonal. 
Basic dynamical properties of the discrete solitons are also discussed in the context of an analytical 
approximation, in terms of an effective Peierls-Nabarro potential in the lattice setting. 
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Introduction. In the past several years, tremendous 
progress has been observed in studies of nonlinear dy- 
namical systems on lattices Q. In a considerable part, 
this development is fueled by the continuous expansion of 
relevant physical applications, including spatial dynamics 
of optical beams in waveguide arrays [2| , temporal evolu- 
tion of Bose-Einstein condensates in deep optical lattices 
0, transformations of the DNA double strand 4], and 
so on. 

A ubiquitous dynamical-lattice system is represented 
by the discrete nonlinear Schrodinger equation [j], 0, Q 
with cubic (x^) nonlinearity. It has been used to model 
a variety of experiments featuring, among others, forma- 
tion of discrete solitons, lattice modulational instability, 
buildup of the Peierls-Nabarro (PN) barrier impeding the 
soliton motion, diffraction management, and soliton in- 
teractions H H H OH EH . 

Substantial activity has also been aimed at lattices 
with quadratic (x^) nonlinearity, boosted, in partic- 
ular, by the recent experimental realization of discrete 
X^ solitons in nonlinear optics A variety of top- 

ics have been studied for such media both theoretically 
and experimentally, including the formation of ID and 
2D solitons pjj, LLJ, LL3 (see also reviews pii Il7j). ob- 
servation of the modulational instability in periodically 
poled lithium niobate waveguide arrays, finite few-site 
lattices [HI, X^ photonic crystals Ull. cavity solitons 
[20j | . and multi-color localized modes |2l|. In addition, 
the same lattice models with the quadratic nonlinearity 
may be used to describe the dynamics of Fermi-resonance 
interface modes in multilayered systems based on organic 
crystals |22j. A variety of solutions have been obtained 
in the latter context [2j . 

A fundamental difference of x^ continua from their 
X^ counterparts [2J is that they feature no collapse in 
2D and 3D cases l25| . which paves the way to create sta- 
ble 2D [2^ and 3D |27j quadratic solitons. On the other 
hand, due to the presence of collapse-type phenomena 
in 2D and 3D x^ continua, lattice solitons in the cor- 
responding discrete setting are subject to quasi- collapse. 



As a consequence, they only exist with a norm ("mass") 
exceeding a certain threshold value [2^, and they are 
strongly localized (on few lattice sites), hence 2D and 3D 
X^ solitons are strongly pinned to the lattice and cannot 
be motile H^. 

The absence of the trend to the catastrophic self- 
compression in the 2D x^ medium suggests that the 
corresponding lattice solitons may be broad and there- 
fore mobile, being loosely bound to the lattice. The aim 
of this work is to investigate the mobility of ID and, 
especially, 2D solitons in x^ lattices. Besides its impor- 
tance for the use of such waves in photonic applications, 
the topic is of fundamental interest by itself, as it will 
reveal a family of mobile solitons in 2D lattices. Thus 
far, the only example of mobility was provided by soli- 
tons in a 2D lattice with saturable nonlinearity 21] (the 
Vinetskii-Kukhtarev model , in which the mobility of 
ID solitons was examined in Ref. HJ). I n this work, 
we identify parametric regions of stable motion of x^ 
solitons on ID and 2D lattices and, for the first time in 
the 2D case, we study anisotropy of the mobility of 2D 
lattice solitons (for propagation off of the principal direc- 
tions of the lattice). First, we will introduce the model 
and develop an analytical approach to solitary waves and 
the respective PN barrier in x^ lattices. Then, system- 
atic numerical results for the soliton mobility in ID and 
2D lattices will be reported. 

The model and analytical results. Following Ref. 
jl5j | . we introduce a system including equations for the 
fundamental- frequency (FF) and second-harmonic (SH) 
waves, il>m, n (t) and m , n (£). In the 2D setting the model 
has the following form (its ID counterpart will also be 
used below): 
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Where A 2 ^ m , n = ^m+l,n + ^m-l,n + ^m,n+l + ^m.n-l 

4w mn is the discrete Laplacian, C\ and C 2 are the FF and 
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SH lattice-coupling constants, and k is the mismatch pa- 
rameter. Equations QJ and © conserve the correspond- 
ing Hamiltonian and the Manley-Rowe (MR) invariant, 

Stationary solutions are looked for as 
{VVnW,</w(£)} = {e~^m,n,e~ 2i ^$m,n}, where 
distributions ^ m ,n, ^m,n are real for fundamental soli- 
tons, and may be complex for more elaborate patterns, 
such as vortices To set discrete solitons in motion, 
one must overcome the above-mentioned PN barrier, i.e., 
the energy difference between static solitons, centered, 
respectively, on a lattice site and between sites. To ob- 
tain an approximate analytical expression for the barrier, 
we consider the continuum limit, in which stationary 
functions \I/ and <1> depend only on the radial variable 
which is the continuum limit of r = y (m 2 + ri 2 ) jC\ 
and obey equations 



(4cj • 



k) $ + C ($" + $'/r) + ^ 2 = 0, 



(3) 



where C = C2/C1, and the prime stands for d/dr. In 
this limit, the soliton is approximated by the following 
ansatz with amplitudes A and 5: 
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where uo and x = (4uo + k) /C must be negative. These 
expressions have the correct 2D asymptotic forms at r^ 

00, {\I>, <£} ~ r _1/2 exp (^-\/R^x} r ) • Substituting the 
ansatz in Eqs. © and demanding its validity at r — > 0, 
we obtain B = (23/3)|cj|, A = (23/3) (4u + k). 

The Hamiltonian corresponding to the axially symmet- 
ric real solutions of the continuum equations is 



H 

-k<$> 2 
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] = 7T / rdr (2cj^ 2 + 4cj$ 2 + 'M' 2 ) , (5) 
Jo 



where the derivatives were eliminated using inte- 
gration by parts and Eqs. ©. To derive 
the PN potential, we apply the lattice discretiza- 
tion to final expression JSJ by defining #i a tt = 
- / / (2u^ 2 + 4cj^> 2 + $^ 2 ) gr(x, y) dxdy, where the 
grid function is 



+00 



+00 



gr(x, y)= Hx-m)S(y-n)= ^ g 2«(^+9y) . 

ra,n= — oo p,q= — 00 

(6) 

In the quasi-continuum approximation (which implies 
small I a; I and |x|), the leading terms in H\ att corre- 
spond to (p, q) = (±1,0) and (0, ±1) and yield the 



PN potential with an exponentially small amplitude, 
U = Uo [cos (27r£) + cos (27777)], where (£,77) are the co- 
ordinates of the soliton's center. An expression for the 
amplitude Uo is simplest in the case of |x| > 2|o;|, 
which corresponds to numerical results presented below 
(with uo = —0.25, x = —0.75); in this case, the sec- 
ond term in ifi a tt dominates. Fitting the slowly vary- 
ing part of the integrand to a Gaussian, we thus ob- 
tain Uo = -a(|cj| 3 /|x|)exp(-37r 2 /(10|x|)), with a = 
(2tt/15)23 2 « 222. 

Numerical Results. In the ID and 2D cases alike, 
we used lattices with periodic boundary conditions, in 
order to allow indefinitely long progressive motion of 
solitons. First, we found standing lattice-soliton solu- 
tions l^m^nj $m?n}j by means of fixed-point iterations. 
Next, dynamical simulations were initialized by applying 
a shove (kick) to those solutions, which corresponds to 
initial conditions 
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where S and 6 determine the size and orientation of the 
shove vector. 

Examples of stable motion and destruction of the ID 
lattice soliton subjected to the shove are displayed in 
Fig. 12 and systematic results, obtained with variation of 
S and C\ = C2, are summarized in Fig. [3 The destruc- 
tion was registered as the outcome if the kicked soliton 
would eventually lose more than 30% of the initial value 
of its MR invariant. For the coupling strength such as 
that corresponding to Fig. [1] there are, practically, only 
two outcomes, mobile waves and wave destruction, ob- 
served with different initial kicks. However, for weaker 
couplings (i.e., stronger discreteness), "localization" is 
also possible: if the shove's strength, 5, is below a lower 
critical value, Scr\ the soliton survives without acquiring 
any velocity. The latter outcome is explained by noting 
that the kinetic energy, E^ n ~ S 2 , imparted to the soli- 
ton by the shove, may be insufficient to overcome the 
PN barrier (2Uq as defined above). Because Uo decays 
exponentially with the increase of the intersite coupling, 
the "localization" region in Fig.[2|is very small. General 
features of the ID situation are: (i) for S < S&\ the 
soliton remains quiescent; (ii) for Sc? < S < 5 cr , the 
soliton sets in the state of persistent propagation; (hi) 
for S > £ cr , the soliton is destroyed. 

We now turn to the 2D setting, which is more inter- 
esting for two reasons. First, as noted above, in the 2D 
case the mobility of lattice solitons is a highly nontriv- 
ial feature, practically impossible in the case of the x^ 
nonlinearity; second, it is interesting to study anisotropy 
of the mobility, i.e., its dependence on the orientation 
of the initial kick relative to the (principal directions of 
the) lattice. Figure [3| shows two examples of stable mov- 
ing regimes: one along the lattice diagonal, and, to our 
knowledge, the first ever example of the motion on the 
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FIG. 1: (Color online) Space-time contour plots of l^m,™! 
and |0 m ,n| 2 for the FF and SH fields (top and bottom pan- 
els) in the ID lattice with periodic boundary conditions, for 
C\ — C2 = 1, uj = —0.25, k = 0.25, and the shove's strength 
S = 0.4 and 3.0 (left and right panels respectively). The 
boosted soliton sets in stable motion in the former case, and 
is destroyed in the latter case. 
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FIG. 2: A diagram in the plane of the coupling strength, 
Ci = C2, and shove's strength, 5, showing different outcomes 
of kicking the quiescent soliton in the ID lattice, for uo — 
—0.25 and k = 0.25. ("localization" means that the soliton 
remains quiescent). 



lattice in an arbitrary direction (neither diagonal, nor 
along the bonds). 

In Fig. 21 we summarize the dependence of the mobility 
properties on the strength, 5, and direction, of the 
initial kick. The kicking of the 2D soliton may result in 
"localization" (no motion at all), in some interval S < 

S£ ] (S < S£ ] « 0.02 in the top left panel of Fig. S|). 
Other generic outcomes again amount to propagation at 
a finite velocity, which depends on 5, and destruction for 
a large supercritical kick S > S cr . 

Particularly noteworthy features, specific to the 2D 
setting, are presented in the top right and bottom left 
panels of Fig. 21 viz., dependences of S cr and velocity in 
the moving regime on 0. These dependences demonstrate 
that the propagation direction along which it is easiest 
to sustain motion (i.e., with larger speeds/range of ini- 
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FIG. 3: (Color online) Same as Fig. Q]and with the same pa- 
rameters, but in the 2D periodic lattice, for the propagation in 
the diagonal (45 degrees, top panels) and off-diagonal (20 de- 
grees relative to the lattice bonds, bottom panels) directions. 
The left panels show trajectories of the soliton's center, while 
the right ones display snapshots of the moving solitons in the 
FF (top) and SH (bottom) fields at t = 0, 30, 50. 
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FIG. 4: (Color online) Features of the soliton motion in the 2D 
periodic lattice, for Ci = C 2 = 1, k = 0.25, u = -0.25. The 
top left panel shows the velocity versus the shove's strength 
S in the propagation along the lattice bonds (at angle = 0); 
the vertical dashed line indicates the value of 5 cr , beyond 
which the soliton retains less than 70% of its initial MR in- 
variant, and is therefore categorized as destroyed. The top 
right panel depicts S C r as a function of the orientation of the 
initial kick, 0. For a given S (S = 0.4), the ensuing velocity 
of the motion is shown versus 6 in the bottom left corner. 
In addition, the bottom right panel shows the analytically 
predicted PN potential. 

tial kicks) on the square lattice is along the diagonal, as 
the motion in this direction can be sustained up to larger 
values of S cr , and is fastest for given S. Both facts may 
be qualitatively explained by the analytically predicted 
PN potential in the following way. Given the nearest- 
neighbor nature of the interactions, in order for the cen- 
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ter of the wave to move along the diagonal, it has to split 
along the two lattice directions and then recombine at 
the site located diagonally across from the initial posi- 
tion. The recurrent small-scale symmetric breakups and 
recombinations (observed in the numerical data) provide 
for an effective propagation along the diagonal direction 
with a minimum PN barrier (see bottom right panel of 
Fig. |3} and thus result in higher propagation speeds. 

We have also examined the situation with C2 < Ci, 
and obtained similar results, but with larger Sc? • In the 
special case of C2 = (no lattice coupling in the SH field), 
we were not able to generate moving solitons, which can 
be easily explained: with C2 = 0, Eq. yields <& m ,n = 
— (4uj + &) _1 (^ m , n ) 2 •> and the substitution of this in Eq. 
JQ makes the model equivalent to that with the cubic 
nonlinearity, where moving 2D discrete solitons do not 
exist. 

Conclusions. In this work, we have examined the mo- 
bility of solitons in ID and 2D lattices with the quadratic 
nonlinearity. We have shown that the solitons feature 
stable motion much more easily than their counterparts 
in ID lattices with the cubic nonlinearity, and they may 
also be mobile on the 2D lattice, where the cubic solitons 



cannot move at all. In the 2D lattice, we have for the first 
time reported a possibility of motion of the soliton in an 
arbitrary direction (neither axial nor diagonal), while we 
have illustrated the interesting differences between prop- 
agating on and off of lattice directions. A qualitative 
explanation for some of these features was provided by 
an analytical approximation for the 2D Peierls-Nabarro 
potential. 

It may be interesting to extend this type of exami- 
nation to other ID and, especially, 2D models, where 
mobile solitons may be expected, such as systems with 
competing nonlinearities (the cubic-quintic model , or 
the Salerno model with competing on-site and inter-site 
cubic terms js3) an d saturable nonlinearity [2^, where 
one may expect significant potential for genuine traveling 
[3lL 134^ . While herein the mobility of solutions for real- 
istic physical purposes in finite lattices was examined, 
the existence of genuinely traveling such solutions is also 
an interesting computational [^H and mathematical j^H 
problem. 
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